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Abstract. We describe an algorithm to compute the local component at p of 
the Coleman-Gross p-adic height pairing on divisors on hyperelliptic curves. 
As the height pairing is given in terms of a Coleman integral, we also provide 
new techniques to evaluate Coleman integrals of meromorphic differentials and 
present our algorithms as implemented in Sage. 



1. Introduction 

For an elliptic curve over Q, the classical Birch and Swinnerton-Dyer (BSD) 
conjecture predicts that a special value of its L-function can be given in terms 
of certain arithmetic invariants of the curve, one of which involves the canonical 
height pairing matrix of a basis of rational points. The p-adic analogue [MT T86] of 
the BSD conjecture makes a similar prediction, with the canonical height pairing 
replaced by a p-adic one |MT83] . These conjectures have natural generalizations to 
abelian varieties. 

The p-adic height pairing was first defined by Schneider [Sch82| for abelian va- 
rieties and was extended to motives by Nekovaf [Nck93 . For Jacobians of curves 
there is a third definition, due to Coleman and Gross [CG89] , relying on Coleman's 
theory of p-adic integration |Col821 ICol851 ICdS88[ lB"es02) . This third definition of 
the height pairing is known to be equivalent to the previous ones |Bes04] . 

For the purpose of numerically verifying p-adic BSD type conjectures, it is im- 
portant to have an effective algorithm for the computation of the p-adic height 
pairing. By the work of Kedlaya [KedOlj and Mazur, Stein, and Tate [MST06] ■ we 
can easily compute p-adic heights on elliptic curves. Our work deals with the next 
logical step, p-adic height pairings on Jacobians of hyperelliptic curves. 

The reason for treating Jacobians of hyperelliptic curves is that we have available 
the Coleman-Gross definition of the height pairing, which is much more concrete 
than previous definitions. The restriction to hyperelliptic curves is made primarily 
so that we may apply the recent algorithm [BBKlOj for the computation of Cole- 
man integrals on such curves, relying in turn on Kedlaya's work on the computation 
of the matrix of Frobenius on hyperelliptic curves [KedOlj . We note that general- 
izations of Kedlaya's work to other types of curves can be applied to generalize the 
results to these curves as well. 

Coleman and Gross give a decomposition of the global p-adic height pairing as a 
sum of local height pairings at each prime. The local heights away from the prime 
p behave in much the same way as local archimedean heights, so the main interest 
lies in the primes above p, where Coleman integration is used. It is this last type 
of local height pairing which we aim to compute. 
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To fix ideas, consider a hyperelliptic curve C over Q p with p a prime of good 
reduction. Then, for D\,D2 <E Div°(C) with disjoint support, the Coleman-Gross 
p-adic height pairing at the prime p is given in terms of the Coleman integral 



for an appropriately constructed differential uj£ )l associated to the divisor D\. This 
last association is not straightforward and relies again on Coleman integration. 

We say next to nothing in this work about the computation of local height 
pairings away from p. As mentioned before, this is not a p-adic problem and 
is shared with the computation of archimedean height pairings. In lBes07| we 
suggested a method for treating this problem, which needed in particular some 
refined estimates of Kausz [Kau99 . In the meantime, we have learned that this 
problem is treated in the recent Ph.D. thesis of Miiller |MullO| . 

The structure of the paper is as follows: in Section [5J we review the work of 
Coleman and Gross [CG89] that constructs p-adic heights for curves, and in par- 
ticular, defines the local contribution at p in terms of a Coleman integral. After a 
brief review of hyperelliptic curves in Section [3] we give, in Sectional an overview 
of Coleman integration, where we discuss some known results for computing these 
integrals and describe a new construction that allows us to compute a broader class 
of integrals that the heights necessitate. We discuss the algorithm for computing 
the local height pairings in Section [5] In Section [BJ we discuss our implementation 
of the algorithm in Sage along with error bounds on our results. We follow this 
in Section [7] with numerical examples illustrating our methods. We conclude in 
Section [5] by posing some questions arising from our work. 
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In this section we review the definition of the Coleman-Gross height pairing. As 
explained in the introduction, there are two required ingredients for making this 
definition: the theory of Coleman integration and a certain choice of a canonical 
form. These will be discussed in detail in a later section. 

Suppose X/ K is a curve defined over a number field K , with good reduction at 
primes above p. To define the height pairing 
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2. The p-ADic height pairing 



h : Div u (X) x Div°(X) Q. 
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one needs the following data: 

• A "global log"- a continuous idele class character 

£ : k*/K x -> Q p . 

• For each v\p a choice of a subspace W v S H^-^X^Ky/Ky) complementary 
to the space of holomorphic forms. 

For the definition of the Coleman- Gross height we must insist that the local char- 
acters £ v induced by £, for v\p, are ramified in the sense that they do not vanish on 
the units in K v . 

From £ one deduces the following data: 

• For any place v \ p we have £ v (0^ ) = for continuity reasons, which 
implies that £ v is completely determined by the number £ v (n v ), where 7r„ 
is any uniformizer in K v . 

• For any place v\p we can decompose 




where t v is a Q p -linear map. Since we assume that £ v is ramified it is then 
possible to extend 

(2.1) log, : -> K v 

in such a way that the diagram remains commutative. 

Let us now describe the height pairing h(D\ , D 2 ) for a pair of degree zero divisors 
D\ and D 2 with disjoint support. The height pairing is a sum of local terms 

h{D l7 D 2 ) = Y / h v (D 1 ,D 2 ) 

V 

over all finite places v. The local terms depend only on the completion at v of K. 
Thus, let A; be a local field of characteristic 0, with valuation ring O, uniformizer 
7r and let F — O/irO be the residue field, with order q. Let C denote the curve X 
over the local field k. We shall assume that C has a fc-rational point and that C 
has good reduction at 7r. 

Let x '■ k* — > Q P be a continuous homomorphism, which is the local component 
of£ 

Proposition 2.1. // char F ^ p, and Di, D2 6 Div°(C) have disjoint support, then 
there is a unique junction (Di,D2) that is continuous, symmetric, bi- additive, taking 
values in Q p , and satisfying 

(2-2) ((f),D 2 )= X (f(D 2 )) 
forfek(C)*. 

Proof. See |CG89l Prop 1.2]. In fact, one has |CG89I (1.3)] 
(2.3) h v (D u D 2 ) =£ v (ir v ) ■ (Di,D 2 ) . 

Here, (Di,D 2 ) denotes intersection multiplicity on a regular model of C over of 
extensions of Di and D 2 to this model. To make this have the required properties 
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one of these extensions has to have zero intersection with all components of the 
special fibre. □ 

We now describe the local contribution at a place v\p. 

Definition 2.2. A meromorphic differential on C over k is said to be of the first 
kind if it is holomorphic, of the second kind if it has residue at every point, and 
of the third kind if it has a simple pole with residue in Z, respectively. 

Recall that the differentials of the second kind, modulo exact differentials, i.e., 
differentials of rational functions, form a finite dimensional fc-vector space of dimen- 
sion 2g. It is canonically isomorphic to the first algebraic de Rham cohomology of 
C/k, H^ R (C/k), which is the hypercohomology group of the de Rham complex 

— > O c — ► Sl c/k — > 

on C. We have a short exact sequence 

(2.4) — > H°(C, Q c/k ) — > HUC/k) — ► H\C, O c ) — > 0, 

where, relying on the description of de Rham cohomology in terms of forms of the 
second kind we have 

• H°(C, f2^ fe ), the space of differentials of the first kind, is identified with 

its image. It has dimension g, and we will denote it H^(C/k). 

• H (C, Oc) also has dimension g and may be canonically identified with the 
tangent space at the origin of the Jacobian of C, J = Pic (C) . 

• H^ R (C/k) has a canonical non-degenerate alternating form given by the 
algebraic cup product pairing 

H l m [C/k) x H\ K {C/k) — ► k 

(MJH) ^ [Mi] U [/u 2 ], 

which can be described by the formula 

(2.5) M U fa] = ^Res P (/i 2 I 

p ^ 

where are differentials of the second kind, with classes and [/12], 

respectively, in H^iC/k) and the sum is over all points in C. The residue 
does not depend on the choice of a particular local integral for fix because 
fj,2 is of the second kind and has no residue at any point. 

We will also need the theory of Coleman integration. Details will be discussed in 
Section [4] but for now it suffices to know that for each meromorphic form uj on C 
and to each D e Div°(C) , the theory allows us to define the integral J D uj € k. In 
particular, for two fc'-rational points P and Q, where k! is any finite extension of 
k, it allows us to define the integral J? uj G k' . This construction is such that as Q 
(or P) varies, this last integral is expressed as a power series in Q which is locally 
convergent and whose differential is uj. For forms that have residues, the Coleman 
integral depends on the choice of a branch of the p-adic logarithm function. We fix 
this choice for the computation of the local height pairing to be the one determined 
in (l2~TTl . 
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Let T(k) denote the subgroup of differentials on C of the third kind. We have a 
residue divisor homomorphism 

Res : T(k) — > Div°(C), Res(w) = ^ Res P w 

p 



where the sum ranges over all closed points of C . That the image is contained in 
Div°(C) is just the residue theorem. By the residue divisor homomorphism, T(k) 
fits into the following exact sequence: 

(2.6) fSl^C/k) ^T(k)—^lBiv a (C) »-0. 

We are interested in a particular subgroup of T(k) whose elements are the log- 
arithmic differentials, i.e., those of the form y for / e k{C)* . We denote this 
subgroup as Tj(fc). Since T t (k) n H^{C/k) = {0} and Res(^) = (/), we deduce 
from the sequence (|2.6p the short exact sequence 



— > H^(C/k) — > T{k)/Ti{k) — ► J(k) — ► 0. 

This sequence has a natural identification with the fc-rational points of an exact 
sequence of commutative algebraic groups over k: 

— > H°(C, nj, /k ) ► E > J > 0, 

where E is the universal extension of J by a vector group and = Gg. Since 

the Lie algebra of -E is canonically isomorphic to Hi R (C/k), the exact sequence 
(|2.4|) is the resulting exact sequence of Lie algebras over k. 

Now as k is p-adic, we will make use of the fact that we have a logarithmic 
homomorphism defined on an open subgroup of the points of any commutative 
p-adic Lie group, G, to the points of its Lie algebra Lie(G). When G — E or 
J, the open subgroup on which the logarithm converges has finite index, so the 
homomorphism can be uniquely extended to the entire group. We denote this 
extension as \og E or logj, respectively. Since the logarithm is functorial and equal 
to the identity on H°(C, 0^ /fc ), we obtain the following: 

Proposition 2.3. There is a canonical homomorphism 

* : T(fc)/lKfc) — > Hl K (C/k) 

which is the identity on differentials of the first kind and makes the following dia- 
gram commute: 

^ H\£(C/k) E(k) *~ J(k) ^ 



OSS. j 



H^(C/k) H\ K {C/k) Hi(C, O c/k ) 0. 

Note that the map *f> takes a differential of the third kind on C to a differential of 
the second kind modulo exact differentials. It can be extended to a linear map from 
the A:- vector space of all differentials on C/k to H^ R (C/k) by writing an arbitrary 
differential v as a linear combination v = Y] ctifii + 7, where //$ is of the third kind, 
on E k, and 7 is of the second kind on C. We then define = a i^il Jl i) + H- 

The definition of the log map "J/ is not very useful for computations. An equiv- 
alent alternative definition has been given in BesOO] . It is based on the notions of 
local and global symbols or indices. 
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Definition 2.4. For u) a meromorphic form and p a form of the second kind, we 
define the global symbol as a sum of local symbols (lu,p)a- We have 

(u,p) = ^2(uj,p)a, 

A 

where 

(2.7) (u, p) A = Res A 

Note that the sum is taken over all points A where either porw has a singularity. 
Each local symbol (uj, p)a is computed in the local coordinates of A, and Z is a 
point that is fixed throughout the calculation of a single global symbol, which sets 
the constant of integration for the first indefinite integral of p. More precisely, we 
note that (|2.7p is computed by rewriting w and p using the local coordinates at A, 
computing the indefinite integral of p (as a formal power series, and in particular, 
not as a Coleman integral) and computing the Coleman integral of p from Z to A. 

The following result |BesOO] reduces the computation of ^ to the computation 
of global symbols. 

Proposition 2.5. Let to be a meromorphic form and p a form of the second kind. 
Then {u},p) = <F(w) U [p]. 

Proof. See |Bes00i Prop 4.10]. □ 

Now recall that we have at our disposal the complementary subspace W — W v . 
It allows us to isolate a canonical form urj with residue divisor D as follows. 

Definition 2.6. For any divisor D of degree on C we let curj be the unique form 
of the third kind satisfying 

res(cjc) = D, *(wd) £ W. 

It is easy to see from the properties of ^ that this indeed uniquely defines the 
form u)d- 

Definition 2.7. The local height pairing is defined by 

h v (D 1 ,D 2 ) := t v (^J uj Di 

(recalling that the supports of D\ and D 2 are disjoint), where t v is the trace map 
determined by the decomposition of £ v . 

Remark 2.8. In certain cases, there is a canonical complement W to H^(C/k) in 
H^ R (C/k). Namely, when C has good ordinary reduction, we may take W to be 
the unit root subspace for the action of Frobenius. 

Some properties of the local height pairing are as follows: 

Proposition 2.9. The local height pairing h v (Di, D 2 ) is continuous and bi-additive. 
It is symmetric if and only if the subspace W of H^(C/k) is isotropic with respect 
to the cup product pairing. Finally, the formula (|2.2[) holds. 

Proof. See |CC89l Prop 5.2]. □ 
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3. HYPERELLIPTIC CURVES 

In this section we review relevant facts about algebraic and computational as- 
pects of the theory of hyperelliptic curves. 

Let us now suppose that C is a hyperelliptic curve. As k has characteristic not 
equal to 2, a hyperelliptic curve C over k of genus g is an algebraic curve given by 
the equation 

(3.1) y 2 = f(x) , f(x) e k[x], 

where / has simple zeros. For our implementation in Sage, we shall further assume 
/ to be monic and deg / = 2g + 1. 

The curve C is singular only at infinity (and non-singular when g = 1). To 
describe the neighborhood of infinity we normalize the curve there and obtain the 
equation at infinity 

(3.2) t 2 = sf ev (s) withx = -, y 



S' " S9+ 1 ' 

where f rcv (s) = s 29+1 f(l/s) is the reversed polynomial. The point at infinity 
corresponds in these coordinates to both s and t being 0. Furthermore, t has a 
simple zero and s has a double zero at this point. 

As is well-known, the first de Rham cohomology of C has a basis consisting of 
the forms of the second kind 

x 1 dx 

oji := — — for i = 0, . . . , 2g - 1. 
We will denote this basis as 

f dx xdx x 2g ~ 1 dx 

(3 - 3) »-{*•-*■■■■■— 

If we make the change of coordinates (|3.2j) we see that these are transformed as 
follows: 

x z dx „_-, jds 

— — i — ^ -s 9 l — . 
2y 2t 

Since s has a double zero at the point at infinity one sees that these forms are 
holomorphic for i = 0, . . . , g — 1 and meromorphic otherwise. 

We finally recall the hyperelliptic involution w defined by w(x,y) = (x, —y). 

4. Coleman integrals 

Here we review the relevant background on Coleman integrals and describe new 
techniques to handle Coleman integrals of meromorphic differentials with poles in 
non-Weierstrass residue discs. This gives us the necessary tools to present our 
algorithm to compute local heights. 

4.1. Differentials of the second kind. The foundational reference for this is 
Col85]. A more expanded version of our presentation can be found in [BB KlOj . 

Let ui be a 1-form, with (w)^ denoting its polar support. For P, Q G C(fc), 
Coleman integration allows us to compute f p 0J £ k. 

Theorem 4.1. Let fi,v be 1- forms on C and P,Q,R G C(k). The (definite) 
Coleman integral has the following properties: 

(1) Linearity: J p (afj, + bv) = a f p (j, + b f p v, for P,Q £ (p)oo U (^)oo- 
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(2) Additivity: Jp fi = Jp fi + J^ [i, for P,Q,R f (m)oo ■ 

(3) Change of variables: If C is another curve and <fi : C —t C a rigid analytic 
map between wide opens then Jp <f>* fj, = J^p } M */ 4 > {P)i <^(Q) ^ {l^)oo- 

(4) Fundamental theorem of calculus: jp df — f{Q) — f(P) for f a meromor- 
phic function f on a wide open subset. 

Proof. See |Col85| Thm 2.3, Prop 2.4, Thm 2.7] for details. □ 

Thus writing 

Ld = df + CqUJo H h C2g-\U)2g-l, 

where / is a function and Ci 6 k, we see that the problem of computing J p lo 

reduces to computing the Coleman integral of a basis differential, Jp u>i. 

For hyperelliptic curves, Coleman integration can be performed numerically. For 
simplicity we describe things only over Q p although things work in complete gen- 
erality (see [BBK10] and also [BeslOj where a slightly different but related version, 
which was implemented in |Gut06j . is described). 

Let us recall the following algorithms (Algorithms 8 and 11, respectively) from 
BBK10J. Throughout, let <j> denote a p-power lift of Frobenius. 

Algorithm 4.2 (Tiny Coleman integrals). 
Input: 

• Points P,Q £ C(Q P ) in the same residue disc (neither equal to the point 
at infinity). 

• A basis differential w,. 

Output: The integral J p Wj. 
The algorithm: 

(1) Construct a linear interpolation from P to Q. For instance, in a non- 
Weierstrass residue disc, we may take 

x(t) = (l-t)x(P)+tx(Q) 

y (t) = vTR^yy, 

where y(t) is expanded as a formal power series in t. 

(2) Formally integrate the power series in t: 

Jp Jp 2y Jo 2y(t) dt 

Algorithm 4.3 (Coleman integration in non-Weierstrass discs). 
Input: 

• The basis differentials {wjj-^g 1 . 

• Points P,Q G C(Q p ) in non-Weierstrass residue discs. 

Output: The integrals {/p Wi}-^ 1 . 
The algorithm: 

(1) Calculate the action of Frobenius on each basis element : 

2g-l 

<j>*{u i )=df i + J2 MijU) r 
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(2) By change of variables (see Remark l4.4j) . we obtain 

Sg" 1 rQ r<t>(P) rQ 

(4.1) V(M-J)«/ Uj = fi(P) - fi(Q) - ut- «h. 

j=0 Jp Jp U{Q) 

As the eigenvalues of the matrix M are algebraic integers of C p -norm p 1 / 2 ^ 
1 (see |Ked01| §2]), the matrix M — I is invertible, and we may solve (|4.1j) 
to obtain the integrals J p uli. 

Remark 4.4. We obtain (|4.1I) as follows. By change of variables, 

r<«Q) 

0*(Wi) 

P 

(4f 4 + X M «tf w i) 

3=0 

-"■'■'./ / ^ /• 

3=0 



MQ) - /i(-P) + >. M« / a." 

P 



Adding + J^q) to both sides of this equation yields 

,0(P) ,Q 2g-l Q 

Ui = / «< + / Wi + - /i(P) + V My / Uj , 

P JP J<t>(Q) j=0 JP 

which is equivalent to (|4.1I) . 

4.2. Meromorphic differentials. The above approach does not work for a mero- 
morphic differential with poles in non-Weierstrass residue discs. We provide a new 
approach for dealing with this case (again we describe things only over Q p but 
things work in general). Let ui be such a differential. As before, if R, S are points 
in the same non-Weierstrass residue disc (different from those discs containing the 
poles of uj), then J s uj is just a tiny integral, which can be computed as in Algo- 
rithm 14.21 Let us now suppose that R, S are in different non-Weierstrass residue 
discs. 

Algorithm 4.5 (Coleman integration: differential with poles in non-Weierstrass 

discs). 

Input: 

• The differential uj with residue divisor (P) — (Q), with non-Weierstrass 
points P,Q <E C(Q P ). 

• Points R,Sd C(Q P ) in distinct non-Weierstrass residue discs. 

Output: The integral f s ui. 
The algorithm: 

(1) Let a — <P*(uj) — pui. Using the methods of Section [5721 compute &(uj). 
Using this, compute *(a) as (f>*(^(uj)) —p^(uj). 

(2) Let f3 be a form with residue divisor (R) — (S). Compute ^f(/3). 

(3) Compute the cup product Vt'(a) U ^(/3) (see Section HTT1 for more details). 

(4) Evaluate the tiny integrals J^gj w and Jr oj- 
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(5) Let S be the set of Weierstrass points and poles of a. For each point A E S, 
expand a J f3 in terms of the local coordinate at A and find the residue at 
A. Compute the sum 

Res^(a f p). 

AGS 

(6) By Coleman reciprocity (see Remark |4.6|) . recover the desired integral: 
/> = T^(*(«)U*(» + X>s. (aft) -J^-C") ■ 



(4.2) 



Remark 4.6. We obtain (14.21) as follows 



rR [-R rR 

(4.3) / a= (/>*lu-p uj 

J S *■* S J- s 



4>(R) rR 
uj — pi UJ 

4>(S) Js 

R ( r S r 4>(R.) 



= (l-p) / W+ / UJ+ w . 

Js \J<t>( s ) J R ) 

As (3 is a form with residue divisor (R) — (S), by Coleman reciprocity |Col89[ 
Theorem 5.2], we compute the LHS of (|4.3I) to be 

Y, Rcs A [P / ^ = *(a) U *(/3) - £ Res A (<* [ fi) , 
A&T \ J / AeS \ J / 

where T is the set of Weierstrass points and poles of j3. We then have 

f » - rb (* <Q > u + J>" ("/")- jC" " C u ) ■ 

Remark 4.7. In practice, the computation in Step[5]of Algorithm 14. 51 is the slowest 
part of the algorithm, as it involves high-precision local calculations over all poles 
of a and all Weierstrass points of the curve. However, since J2 TeU ResT(a) = 
in each residue disc U, for the Weierstrass discs, we do not need a constant of 
integration. For the poles of a, we may choose one constant of integration within 
each residue disc. More precisely, if P and Q are in separate residue discs, we 
compute 

trfe(x(P!))/fc / P j 



Y Res A J P^j =tok{x{P 1 ))/k ^Rcs Pl J pj^j ^ 

Y R- e SA (a J ' pj =tr k{x{Ql))/k \Res Ql (a J Pjj = tr k{x{Ql) y k (- J (i\ , 



where P\ (resp, Qi) is a root of a in the residue disc of P (resp Q). 
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5. The local height pairing at primes above p 

In this section we will explain the algorithm that computes the local height 
pairing at a prime above p for degree zero divisors on the hyperelliptic curve C . 
Recall that we have as additional data the complementary subspace W and the 
character \ from which we deduce a branch of the logarithm to be used in Coleman 
integration and the trace map t v (we keep the subscript v at some places for clarity, 
even though it now serves no purpose). 

Let Di and D2 be two divisors of degree on C. Our main algorithm computes 
the local height pairing h v {D\,D2). It may be described in two steps 

• Compute the height pairing in the case where D\ and D2 are anti-symmetric 
with respect to the hyperelliptic involution (Algorithm 15.81) 

• Compute the height pairing in the general case using the first case (Algo- 
rithm 15.71) . 

Before discussing either algorithm, we begin with some general notes about the 
representation of divisors on hyperelliptic curves (see |Kob981 App §5-6]). 
Recall that a divisor of degree on C may be written in the form 

D = ^ rriiPi - m^j (00). 

Definition 5.1. A divisor D as above is called semi-reduced if the following con- 
ditions are satisfied: 

• rrii > 

• If Pi is in the support of D, then — P,; is not, unless Pi = —Pi in which case 
mi = 1. 

A semi-reduced divisor is called reduced if in addition 

One may represent a semi-reduced divisor D by a pair of polynomials a(x), b(x) 
with deg(6) < deg(a) such that 

• The projection of ^ rriiPi on P 1 is the zero divisor of a(x). 

• b(x) is an interpolation polynomial with the property that for Pi = (xi,yi) 
we have fr(xj) = yi- 

The condition that D is reduced is equivalent to having deg(a) < g. 

Remark 5.2. One can associate b uniquely to the divisor by insisting that a(x)\(b(x) 2 — 
f{x)). This would be less important for us and there are cases we may not achieve 
this. 

Definition 5.3. Let us denote by (a, b) the semi-reduced divisor determined by 
the pair of polynomials a and b and call (a, b) the standard representation of the 
divisor. 

It is known that any degree zero divisor D on C is equivalent to a unique reduced 
divisor. Furthermore, the reduction is effective. More precisely, passing from an 
arbitrary divisor to a semi-reduced one is just a question of adding or subtracting 
divisors of functions pulled back from P 1 while passing from a semi-reduced divisor 
to a reduced divisor has an effective algorithm described in |Kob981 App, Alg 2 
and Thm 7.2 ]. Since we know our height pairing satisfies (12.21) by Proposition ^. 9[ 
which easily allows to pass from a divisor to an equivalent divisor in the pairing 
there is no harm in assuming that our divisors are reduced. 
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Unfortunately for us, reduced divisors are not sufficient. The reason is that since 
they always have a component at infinity, two such divisors cannot have disjoint 
support unless one of them is trivial. For this reason we will work with the difference 
of two reduced divisors. 

Definition 5.4. The divisor denoted (a, 6) — (c, d), where a and c are polynomials 
of the same degree < g, stands for the difference of the reduced divisors defined by 
(a, 6) and (c, d). 

We always assume that the two divisors defined by (a, b) and (c, d) have no 
common components. If there are common components they can be cancelled out. 

We will mostly work with antisymmetric divisors. Given any zero divisor D, 
the divisor D — w*D is antisymmetric. Conversely, any antisymmetric divisor is 
obtained in this way. It follows easily that any antisymmetric divisor is equivalent 
to D — w*D for a reduced divisor. There may be several representations in this 
form for a given divisor, however, there is just one containing no points Pj with 
Pi = —Pi- In the representation (a, b) for D this is equivalent to having a prime to 
/■ 

Definition 5.5. The standard representation of an anti-symmetric divisor is 

[a, b] := D - w*D 

with D a reduced divisor given in the form (a,b). It is a reduced standard repre- 
sentation if a is prime to /. 

Note that with anti-symmetric divisors the standard representation suffices be- 
cause they do not have components at infinity. Thus, [ai,&i] and [02,62] have 
disjoint support if and only if a\ and 02 are relatively prime. 

We now describe Algorithm 15.71 For any divisor D we have a decomposition, 
with rational coefficients 

(5.1) D = ^D+ + ^D-, D + = D + w*(D) , D~ = D — w*(D). 

Lemma 5.6. Suppose D is given by the representation (a, 6) — (c, d). Then, the 
divisor D~ is just [a ■ c, e] in terms of Definition \5.b\ where e is obtained by solving 
a Chinese Remainder Theorem problem to be congruent to b modulo a and to —d 
modulo c. 

Proof. If a is prime to c this is clear. In general, suppose {x — a) has multiplicity 
m in a and n in c. We may assume m > n. Since we are assuming (a, b) and (c, d) 
have no common components, it follows that D — w*{D) is going to have the two 
summands (m + n)[(a, 6(a)) — (a, —6(a)], and to + n is indeed the multiplicity of 
ac in a (it seems though that we can only solve the Chinese remainder problem 
modulo the least common multiple of a and c). □ 

On the other hand, D + is nothing but the divisor of the rational function 2j^t 

c \ x ) 

considered as a rational function on C. It follows from the fact that (|2.2[) is satisfied 
by Proposition 12.91 that for any E 6 Div°(C) we have 

(5.2) K{D + ,E) = h v (E, D+) = X fyE)) , 

where - (E) means as usual the product of the values of - on the x-coordinates 
of the points making up E with the appropriate multiplicities. An easy consequence 
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of this formula is that 

(5.3) h v (D+,E) = h v (E,D+) = if E is antisymmetric. 

Consider now two divisors Di and D 2 in Div (C). Decomposing into plus and 
minus parts it follows from (|5.3I) that 

(5.4) K(D U D 2 ) = Ui v (D+,D+) + Ui v (D^,D 2 ). 

The first term can be computed using (|5.2[) , while the second term is a height pairing 
between anti-symmetric divisors. This immediately gives the following algorithm. 

Algorithm 5.7 (p-adic height pairing for general divisors). 
Input: 

• The subspace W, branch of logarithm and trace map t. 

• Divisors D\ and D 2 with disjoint support given as 

Di = (ai,6i) - (ci,di) 
D 2 = {a 2 ,b 2 ) - (c 2 ,d 2 ). 

Output: The local height pairing h v (Dx, D 2 ). 
The algorithm: 

(1) Compute expressions for the divisors Dj" and D 2 using Lemma T5. 61 

(2) Compute using ()5.2[) . 

h v (D+,D+) = X ((a 1 /c 1 )(D+)). 

(3) Compute, using Algorithm 15.81 the local height pairing for anti-symmetric 
divisors h v (D± , D 2 ). 

(4) Substitute in (|5.4|) to obtain h v (D 1 ,D 2 ). 

We next turn to Algorithm l 5 ■ 8 1 for the case of anti-symmetric divisors D\ and D 2 . 
First of all, we have to introduce yet another decomposition. The algorithm behaves 
differently with respect to parts that reduce to the Weierstrass points (Weierstrass 
divisors) and those which do not. We can decompose a divisor D into the sum of 
its Weierstrass part D w and its non- Weierstrass part D nw . Then, in a similar way 
to (|5.4p we have the decomposition 

(5.5) h v (D 1 ,D 2 ) = h v (A , DJ) + h v (DY, Df) + K (D™, D™). 

We now give the algorithm and later we discuss each step with some more details. 

Algorithm 5.8 (p-adic height pairing for anti-symmetric divisors). 
Input: 

• The subspace W , branch of logarithm and trace map t. 

• Anti-symmetric divisors D\ and D 2 given in standard representation [a\, 6i), 
{a2,b 2 ). 

Output: The local height pairing h v (Di, D 2 ). 
The algorithm: 

(1) Compute the cup product matrix for a basis of H^ R (C / k) , as in Sect ion [5~T1 

(2) Compute Df and Df w . 

(3) Write down forms and with residue divisors and -D" w , respec- 
tively. 
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(4) Compute the form ojdy an d a holomorphic form rj such that 

uo Dr = v ™ - v . 

(5) Compute the tiny Coleman integral h v {D\, D™) = t(J D „(ujrJY + i>f w — rj)). 

(6) Compute the Coleman integral h v (DY , -D" w ) = ^(Jd™ w ^Dfl- 

(7) Compute the Coleman integral h v (Df w , Df) = t(f Dim ^ nw - J D „ W rj). 

(8) Compute /i„(Di,J5 2 ) using (|5.5p . 

We now add some further details on each step in this algorithm. 

5.1. Computing cup products. We first compute the cup product between any 
two elements of the standard basis (|3.3|) for H^ R (X/k). This is easily done using 
the formula (T2~5j) . 

We can be a bit more precise as follows. 

Definition 5.9. The cup product matrix associated to C with respect to B is the 
2gx2g matrix with entry <nj given by the cup product of differentials [a;j_i]U[caj_i], 
normalized so that U [wj-i] = Res(u;j_i J uii-i) . 

By computing in the local coordinates at infinity, we may record the following 

Lemma 5.10. The cup product matrix for C with respect to B satisfies the following 
properties: 

(1) Anti-diagonal elements are given by the sequence 

111 1 1 



2g - 1 ' 2g - 3 ' 3 3 2g - 1 

(2) Entries above the anti-diagonal are 0. 

(3) Diagonal elements are 0. 

Example 5.11. The cup product matrix with respect to B for an elliptic curve is 

^ ) ' ^°^ e su t> s P ace spanned by is isotropic. In particular, for 

genus 1, we may take H^o{C /k) spanned by luq and W spanned by and we need 
not require p to be a prime of ordinary reduction for the pairing to be symmetric. 

5.2. The map \t. We compute ^ of a differential ui by writing 

= C CJ H h C 2s _lCJ2 9 -l, 

and solving for the coefficients c,-. This is done by considering a linear system 
involving global symbols and cup products: 

2g-l 

Recall that as in Definition 12.41 we calculate the global symbol as a sum of 
local symbols, each of which involves a Coleman integral and a calculation in local 
coordinates. This computation is actually much simpler: 

Proposition 5.12. Suppose uj is a form of the third kind with residue divisor D 
which does not contain oo. Then we have 



{uj,uji) = uji+ Res c 
Jd 

where the residue at oo is computed by taking any antiderivative of uji 
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Proof. The sum of local symbols is over all points where either ui or uii has a 
singularity. These are the points in the support of D and possibly the point oo. 
Since ui has a simple pole at each point P in the support of D, the local symbol is 
simply the multiplicity of D at P times uji (where Z is a fixed point throughout 
the global symbol calculation). Summing over all points gives J D ^>i- On the other 
hand, since we are assuming that uj is holomorphic at oo, the choice of the constant 
of integration for uji at oo does not matter. □ 

Now letting N denote the cup product matrix, we have 

/ (uj,uj ) > 
= iV" 1 : 

5.3. Decomposing a divisor D into D w and D nvr . This is very easy to do. 
When the divisor is given in standard form [a, b] one just reduces a modulo the 
prime ir, picks up the part that reduces to Weierstrass points by taking the greatest 
common divisor with the reduction of /, and then applies a Hensel lift to get the 
factor a w of a corresponding to the points reducing to Weierstrass. Then we have 
a nw = a/a w and the divisor decomposition is deduced from that. 

5.4. A form with the required residue divisor. This is an easy task with the 
following: 

Proposition 5.13. Let the anti- symmetric divisor D be given in standard repre- 
sentation [a, b]. The differential form 

a' (x)b(x)dx 
a(x)y 

has simple poles and its residue divisor is D. 

Proof. Suppose that (a, b) = J2 rn iPi~{J2 rn i) (°°)- We can writer = (b/y) dlog(a). 
The form dlog(a) has simple poles at ±Pj with residue rrii while b/y has value 1 at 
Pi and value —1 at — Pj. On the other hand we can also write uj = (a'6/a)(|^) and 
since || is holomorphic, it follows that uj has no poles where a does not vanish. 
Finally, it is easy to see that ui does not have a pole at infinity. □ 

5.5. Finding ojd for a Weierstrass divisor D. Suppose we have already writ- 
ten down a form whose residue divisor is D. Since the singularities of the form 
are contained in the Weierstrass residue discs, it is amenable to the reduction al- 
gorithm done in Kedlaya's algorithm [KedQlj. This means that we may compute a 
representation of u as a linear combination of the basis B plus an exact differential 
dg. Since it follows that ty(uj) is just the above combination of basis elements, we 
need only subtract the appropriate combination of holomorphic basis elements to 
make it reside in W . 

5.6. Finding uio for a non- Weierstrass divisor D. We start with a form ui 
whose residue divisor is D. We compute 9(oj) as in Section I5T21 All that is left to 
do is to let rj be the projection of ty(oj) on along W. 
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5.7. Integration when D2 reduces to Weierstrass points. Suppose that the 
points of D2 reduce to the Weierstrass points. Since D2 is anti-symmetric, this 
means that in computing the integral we need to take the sum of differences over 
pairs of points ±P which reduce to the same point. This is a sum of tiny integrals. 

5.8. Computation when D\ is Weierstrass and D2 is not. In this case, we 
are given, by the previous reduction, the form ujd 1 as a combination cti^i + dg. 
Thus f D ujd 1 = J2 a i Jd 2 w i + 9(^2). Since the points of D2 are in the domain 
where the integrals of the Ui may be computed, this is a standard computation. 

5.9. Computation when both divisors are non- Weierstrass. In this case, 
lud 1 is given in the form v\ — rj, where v\ is a form with residue divisor D\ and rj 
is holomorphic. The integral of 77 poses no problems while that of v\ is discussed 
in Subsection 14.21 

6. Implementation notes 



In this section, we discuss the choices made in our Sage [S + 10 implementation 



and give error estimates on the precision of our results. We work over Q p with a 
precision of n digits; note that if one desires an answer with n digits of precision, 
one has to start with a larger working precision, as seen below. 

We only discuss the computation for anti-symmetric divisors, as the extension 
to general divisors is trivial, as discussed in Section [5] Furthermore, our imple- 
mentation assumes that the divisors are of the form (P) — (-P) for a Q^-rational 
point P, as it is then quite easy to consider cases when the divisor is a sum of such 
expressions. Finally, all computations are done with respect to a particular choice 
of the complementary subspace W, which we describe below. 

6.1. Auxiliary choices. Our algorithm relies on the splitting 

Hl K {C/k) ~ H$£(C/k) © W, 

which allows us to write 

log(fi) = r) + log(u Dl ), 
where rj is holomorphic, and log(w£>i) G W. 

As noted in Example l5.11[ when g = 1, we may simply take H\ R {C/k) spanned 
by loq and W spanned by u>i. Note that for genus g > 1 it does not suffice to 
take W spanned by co g , Wg+i, ■ ■ ■ , ui2g-i, as the resulting subspace is not necessarily 
isotropic. While the local height is independent of basis, it is not independent of 
the choice of W. For the local height to be symmetric, it is necessary that W be an 
isotropic subspace. For g > 1, we thus further require that p be a prime of ordinary 
reduction, so that W can be chosen to be the unit root subspace. Generalizing 
MST06] . we compute a basis for W as follows: 

Proposition 6.1. Let n be the working precision in the underlying base ring Z p; 
so that all computations are done modulo p n . Let Frob denote the matrix of a 
p-power lift of Frobenius, as acting on the standard basis B of H^ R (C/k). Th 
{Frob" u g , Frob" io g+ i, . . . , Frob" W2 9 -i} is a basis for W. 



en 



Proof. With the integral structure provided by crystalline cohomology, it is well 
known that Frob maps the holomorphic forms to p times the integral structure. 
Thus, with W the unit root part decomposing a vector v into uj + rj with u> e W 
and rj holomorphic it is easy to see that Frob" v is in W + p n times the integral 



COMPUTING LOCAL p-ADIC HEIGHT PAIRINGS ON HYPERELLIPTIC CURVES 17 



structure. In other words, up to the prescribed precision, Frob™ v lies in W. On 
the other hand, Frob is invertible so starting with g independent vectors modulo 
the holomorphic differentials one gets g independent vectors in W. □ 

6.2. Precision. Broadly speaking, the p-adic precision of a local height depends 
on two types of calculation: 

(1) Coleman integrals of basis differentials (or otherwise "nice" differentials - 
e.g., holomorphic in the discs corresponding to the limits of integration) 
and 

(2) expansion of local coordinates at a point . 

Each key step of the algorithm in Section 14.51 can be categorized as depending on 
one or both of these: 

• ^>(lo) needs the cup product matrix (local coordinates) and local symbols 
(Coleman integrals of basis differentials) 

• J 77 is a sum of Coleman integrals of basis differentials 

• J lo is defined in terms of 

— tiny integrals (Coleman integrals of a "nice" differential), 

— sums of residues of Laurent series (local coordinates) , and 

— *(a), (as above). 

6.2.1. Precision of Coleman integrals. For more details on the rigorous computation 
of Coleman integrals, see |BBK10j . We recall the following results ( |BBK101 Props 
18-19]): 

Proposition 6.2. Let Jp lo be a tiny integral in a non-Weierstrass residue disc, 
with P, Q defined over an unramified extension of Q p and accurate to n digits of 
precision. Let (x(t),y(t)) be the local interpolation between P and Q defined by 

x(t) = x(P)(l -t) + x(Q)t = x(P) + t(x(Q) - x{P)) 

y(t) = ^JiMt)). 

Let lo = g{x, y)dx be a differential of the second kind such that h(t) = g(x(t), y{t)) 
belongs to 0[[t]\. If we truncate h(t) modulo t m , then the computed value of the 
integral Jp lo will be correct to min{n,m+ 1 — [log p (m + 1)J} digits of (absolute) 
precision. 

Proposition 6.3. Let Jp lo be a Coleman integral, with lo a differential of the 
second kind and with P,Q in non-Weierstrass residue discs, defined over an un- 
ramified extension of Q p , and accurate to n digits of precision. Let Frob be the 
matrix of the action of Frobenius on the basis differentials. Set B — Frob* — /, and 
letm = v p (det(B)). Suppose the relevant tiny integrals have series expansions trun- 
cated modulo t n ~ 1 . Then the computed value of the integral J p lo will be accurate 
to n — max{m, [log p nJ} digits of precision. 

We will now carefully review the precision of each of the objects we computed, 
as an expansion of the overview in Section 16.21 Let w be a differential with residue 
divisor 

D i = (P) (-P) 
and /3 a differential with residue divisor 

D 2 = (Q) - (-Q). 
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The precision of ^(uj) (and ^f(ft)) just depends on the Coleman integral involved, 
as the residue can just be read off of the differential. 

After computing ty(uj) with respect to the standard basis of H^ R (C), we fix a 
splitting of H\ K {C/k) into H^(C/k) © W, which gives rj and ujd 1 = ui — rj. Since 
the height pairing is given by J u>d x , we need to compute the integrals J oj and J r\. 

The integral J g is just a linear combination of the integrals of holomorphic basis 
differentials. On the other hand, the integral of uj requires the computation of ^(ct), 
^-R- es ( a I P)i an d the tiny integral Jq ^ lj. As before, the tiny integral is 
computed with precision as above. 

Since a — (jfoj — puj, we may write ^(a) in terms of things we have already 
computed, namely ^(a) = Frob(\D r (c<;)) — p^f(uj). So need not do more work here. 
However, the precision of ^ Res(a J ft) merits further discussion, as we must con- 
sider its representation in local coordinates. 

6.2.2. Precision of local coordinates. Computing with local coordinates is crucial 
to the algorithm. More precisely, for any point P, we must construct power series 
x(t),y(t) for a local parameter t such that P = (x(0),y(0)). To explicitly compute 
with power series, we need to know where (t-adically) it is acceptable to truncate 
them. 

6.2.3. Precision: Cup product matrix. The first instance this problem arises is in 
the computation of the cup product matrix. Since Uj(cJj) = 2(g — j) — 2, which is 
minimal for j = 2g — 1, we have minu t (wj) = — 2g. Thus it suffices to compute 
each basis differential oJk to a precision of t 2g . Consequently, we compute x(t) to a 
precision of £ 2 ( 2 f _1 ) and y(t) to a precision of t 29 " 1 . 

6.2.4. Precision: ^ Res(a J /?). Now we consider ^^Res^a J /3), where the sum 
is taken over all points A in the set {poles of a, Weierstrass points of C}. We begin 
by looking at the expansion for a — 4>*lj — pus. 

If A is a non- Weierstrass point defined over Q p , then computing a in local 
coordinates is unnecessary, as a just has a simple pole at A. So Resect / P) = 
ResA ot, and we may simply read off the residue: 

• A — P Teichmuller: 1—p 

• A = P not Teichmuller: —p 

• A — —P Teichmuller: p — 1 

• A = — P not Tcichmiillcr: p. 

In the case where A is non- Weierstrass and defined over an extension of Q p , we 
have to compute Coleman integrals and local coordinates, so we must study the 
precision of both, as given by the following corollary of Proposition 16.21 

Corollary 6.4. Let A be non- Weierstrass, defined over degree p extension k' of 
Q p . Let Ua denote the residue disc of A, and let B be a non- Weierstrass point in 
Ua defined over Q p . Suppose a working precision of n p-adic digits (so that A has 
precision pn in a uniformizer n). Let ft be written in terms of the local coordinate 
(x(t),y(t)) at B, so that ft — h(t)dt with h(ti) truncated modulo t mp . Then the 
residue X^5g(7 A Res(o: J ft) has min{n,pm + 1 — |_log„(pm+ 1)J} digits of precision. 

Proof. If A is non- Weierstrass and defined over an extension K of Q p , then the 
contribution of a to the residue calculation just depends on the disc of A: +1 if A 
is in Up, and —1 if A is in U-p, with the constant of integration fixed in each residue 
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disc. However, in this case, since we have multiple poles in each residue disc, we 
must compute Res(a(/ z /3), where the integral of (5 is definite, taking into account 
a constant of integration chosen for each residue disc. More precisely: suppose we 
are working in the residue disc of P, and say A is defined over a degree p extension 
of Q p . Note that we must compute the local coordinates (x(t),y(t)) at P with a 
precision of at least t pm . As the interpolation from P to A is linear, we merely 
make a linear substitution 

x(t) := x((x(A) - x(P))t) 
y(t) :=y{{x{A)-x{P))t). 

This new x(t),y(t)) is used to compute the tiny integral of /3 from P to A, the 
result of which has precision min{n,pm + 1 — [log p (pm + 1)J}. Taking the trace 
from K to Q p accounts for the other poles of a in the disc of P. □ 



Finally, in the case where A is a finite Weierstrass point, we have to compute in 
the local coordinates of A. (Note that we need not compute the residue at (0, 0) if 
on the curve or at infinity.) 

Proposition 6.5. Let a be above and let A be a finite Weierstrass point not equal 
to (0,0). Let (x(t),y(t)) represent the local coordinates at A. Then to compute 
Res(a / /3) with n digits of p-adic precision, we compute (x(t),y(t)) to t 2pn ~ p ~ 3 . 

Proof. We have 

a = 4>*lj — pu 

y(P)px p ~ 1 dx py(P)dx 
= </>(y)(xP - x(P)) ~ y{x-x{P)Y 

where 

J_ = _ p ^/-i/2\ (f(xp)-f(xyy 
m y ££\ * J f{*) pi 

For Res(a / /3) to have n digits of p-adic precision, we must compute n terms of 
the binomial expansion of . 

Recall that for a finite Weierstrass point (a, 0), we have 



x(t) =a+ -?—t 2 + 0(t 4 ) 
9\ a ) 

y(t) = t, 
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where g(x) = Note that by hypothesis, a ^ 0. We compute the i-adic 

valuation of a: 

Vt(a) = vt(4>*ui) since u> only contributes higher-order terms 
/ y(P)px p ~ 1 dx \ 



Vl. 



U(y)(xP - x(P)) J 
1 



l + v t 



My) 

1 -pv t (y) + (n - l)v t 



( X p ± x{P)) 

fix?) - f(x)P 



l-p + (n-l) 



f{x)P 

(2-2p), if« t (/(sf)-/(o;)P)>0 
— 2p, else 



Thus we have 

v t (a) 



p-2pn + 2n-l, if v t (f(x p ) - f(x)P) > 
p — 2pn + 1, else. 



As p — 2pn + 2n — 1 > p — 2pn + 1 for n > 1, we have Vt(a) > p — 2pn + 1. Set 
m = 2pn — p — 1. Since we want Res(a J /3), we need Vt(a J 0) > —1, so we must 
compute /3 to at least t m ~ 2 . To get this precision, we must in turn compute with 
x(t),y(t) to this precision. □ 

7. Examples 

Here we provide some examples of our algorithms. 

7.1. Local heights: genus 2, general divisors. Let C be the genus 2 hyperel- 
liptic curve 

y 2 = - 23a; 3 + 18a; 2 + 40a; = (a; - 4)(x - 2)x{x + l)(x + 5) 
over Qn, and let 

D i = (P) - (Q) 
D 2 = (R) - (S), 

where P = (-4,24),Q = (1,6), R= (5,30),S= (-2,12). We describe how to use 
Algorithm 15.71 to compute the local contribution at p = 11. 
We see that 

'x - x(P) 



D+ = div 



x - x(Q) 
x — x(R) 
v x — x(S) 
Furthermore, we have 



D- 



[(P)-(-P)} + K-Q)-(Q)l 



DJ = div 



\hu(D+, D+) = \log 



D^ = [(R)-(-R)] + [(-S)-(S)}. 



x-x(P) fJ ^_ 



x - x(Q) 



1 



log 



x{R) -x(P)\ (x{S)-x{P) 



x(R)-x(Q) J \x(S)-x{Q) 



2 • 11 + 9 • ll 2 + 7 ■ ll 3 + 2 • ll 4 + 0(11 5 ). 
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Now we compute, using Algorithm I5.8[ the contributions from antisymmetric 
heights (details of which are in Subsection 17. 2[) : 

hu((P) - i-P), {S) - (S)) = 9 • ll- 1 + 5 + 6-11 

hn{(P) - {-P), (R) - (-R)) = 6 • II" 1 + 10 + 7 • 11 + 6 ■ ll 2 + 3 • ll 3 + 7 • ll 4 

M(-Q) - (Q). ( R ) - (~ R )) = 8 • 11-1 + 



8-ll 2 + 9-ll 3 + 3-ll 4 + 0(ll 5 ) 

0(11 5 ) 
0(11 5 ) 



5 + 7 • 11 + 10 • ll 2 + 3 • lr + 7 • ll 4 
M(-Q) ~ (-O)i (~ S ) - ( S )) = + 8 + 7 • 11 + 2 • ll 2 + 7 • ll 3 + 8 • ll 4 + <3(11 5 ), 
which gives 

h n (D^D^) = ^(M(P) - (--P). (-#) - 

^n((P)-(-P),(i?)-(-i?)) + 

fcn((-Q) - (Q), (R) - (-R))+ 
h u ((-Q) - (-QU-S) - (S))) 



Finally, we have 



6- ll~ 1 + 7 + 4- 11 + 4- ll 2 + 3 • ll 3 



\h u {Dl, Dt) + \h lx {D^,D 2 



11 4 + 0(11 5 ). 



= 6-11" 



7 + 6- 11 + 2- ll 2 + 4- ll 4 



0(11 5 ). 

7.2. Local heights: genus 2, antisymmetric divisors. Keeping notation as in 
Subsection 17.11 we describe in more detail how to use Algorithm 15.81 to compute 
the local contribution for one of the antisymmetric divisors: 

h 11 ((P)-(-P),(R)-(-R))- 

For ease of notation, let us call these divisors 

Dp = (P) - (-P), D R = (R) - (-R). 
With respect to the standard basis B, the cup product matrix is 

\ 



N = 














1 

3 










1 










-1 





23 
3 


V 


1 

3 





23 
3 






/ 



Let vp be a differential with residue divisor Dp: we have vp 



2idx 
V(x+4) ■ 



We compute $f(vp) with respect to the basis {wo, uji, Frob" u>2, Frob™ W3}: 
( 



ll- 1 + 9 + 6 • 11 + 3 • ll 2 + 7 • ll 3 + ll 4 + 0(11 5 ) \ 
7 • ll" 1 + 1 + 4 • ll 2 + 2 • ll 3 + 8 • ll 4 + 0(11 5 ) 
7 + 9 • 11 + 7 • ll 2 + 4 • ll 4 + 8 • ll 5 + 0(11 6 ) 
\ 2 + 2 • 11 + 8 • ll 2 + 6 • ll 3 + 7 • ll 4 + 2 • ll 5 + 0(11 6 ) J 

This implies that 

771 = 5 • ll -1 + 6 + 3 • 11 + ll 3 + 7 • ll 4 + 0(11 5 ). 



To integrate i/p, we compute several quantities. Noting that a 
and that the W map is Frobenius equivariant, we have 

fr(a) = ^{(j)*Vp -pvp) = 0*(*(i/p)) -p^(pp). 



' Vp — pvp 
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In particular, this makes the computation of <j>* (^(vp)) rather easy, as we have 
already computed ^(v P ), and all that is left to do is multiply by the matrix of 
Frobenius. We find that 

/ 6- 11 + 5- ll 2 + 2- ll 4 + 9- 11 5 + 0(11 6 ) \ 
. . _ 2 • 11 + 10 • ll 2 + 8 • ll 3 + 6 • ll 4 + 2 • ll 5 + 0(11 6 ) 
V W " 4. H + 6- ll 2 + 2 • ll 3 + ll 4 + 9 ■ 11 5 + 0(11 6 ) 

\ 3-ll + 2-ll 2 + 8-ll 3 + 2-ll 4 + 4-ll 5 + 0(ll 6 ) / 

We wish (5 to have residue divisor Dp, so let j3 — ^jz^y- Then 

U *(/3) = 6 + ll 2 + 9 • ll 4 + 5 • ll 5 + 0(11 6 ). 

To compute ^Res(a J /?), we must sum over all Weierstrass point and poles of 
a. Recall that within a single residue disc, J^Res^a) = 0. Now computing the 
action of ^ on this differential is slightly more complicated, since instead of just two 
non- Weierstrass poles, we have 2p = 2-11 non- Weierstrass poles: those points in 
the residue discs of P and — P with x-coordinate Cn( — 4) 1 / 11 (where j = 0, . . . , 10). 
This means we must work over the splitting field L_4 = Qn(Cn, (— 4) 1//n ) of 
x 11 +4 over Qn to compute the local symbols. Since each set of pth roots is Galois 
conjugate, working over L_4 yields 

^2{v 1 ,u} i ) Pj = tr L _ 4/Qll ((^i,w J ) Pl ), 

3 

where Pj is the point in the residue disc of P with x-coordinate Cu( — 4) 1 / 11 . We 
have the following contribution from the disc of P: 

10 • 11 + 9 • ll 2 + 4 • ll 3 + 3 • ll 4 + 4 • ll 5 + 0(11 6 ), 
and the total contribution from non- Weierstrass points is twice this, or 

9 • 11 + 8 • ll 2 + 9 • ll 3 + 6 • ll 4 + 8 • ll 5 + 0(11 6 ). 
Meanwhile, the sum of contributions from all Weierstrass discs is the following: 
11 + 4- ll 3 + 6- 11 4 + 11 5 + 0(11 6 ). 

We compute the tiny integral 

4>{R) 

v x = 8 • 11 + ll 2 + 8 • ll 3 + 2 • ll 5 + (9(11 6 ). 

£ 

Putting all of this together, we have 

hu{D Pl D R ) = 6- ll" 1 + 10 + 7- 11 + 6 • ll 2 + 3 • ll 3 + 7 • ll 4 + 0(11 5 ). 
As a consistency check, we also compute hu(Dp, Dp). Here we have 

/ ^ = 2 + ll 3 + 10-ll 4 + 4-ll 5 + O(ll 6 ) 

J Dp 

and 

/ rj R = 5 • ll" 1 + 2 + 3 • 11 + 4 • ll 2 + 8 • ll 3 + 2 • ll 4 + 0(11 5 ), 

J Dp 

which gives 

hu{D R , Dp) = 6 • ll" 1 + 10 + 7 • 11 + 6 • ll 2 + 3 • ll 3 + 7 • ll 4 + 0(11 5 ), 
illustrating symmetry of the local height pairing. 
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7.3. Global heights: genus 1. We give an example of our implementation in 
genus 1, which allows for comparison of global heights via the algorithm of Mazur- 
Stein-Tate. 

Let C be the elliptic curve 

y 2 = x 3 — 5x, 
with Q = (-1, 2), R = (5, 10), so that 

(Q)-(-Q) = (#)-(-#)= = : p. 

We compute the 13-adic height of P: 

• Above 13, the local height h 13 ((Q) - (-Q), (R) - (-R)) is 

2 • 13 + 6 • 13 2 + 13 3 + 5 • 13 4 + 0(13 5 ). 

• Away from 13, the only nontrivial contribution is at 3, which is 2 log 3 (by 
work of Miiller). 

• So the global 13-adic height is 12 • 13 + 4 • 13 2 + 10 • 13 3 + 9 ■ 13 4 + 0(13 5 ). 

We compare this to Harvey's implementation |Har08) of the Mazur-Stein-Tate al- 
gorithm in Sage: 

sage: C = EllipticCurve( [-5,0] ) 

sage: f = C.padic_height(13) 

sage: f (C(9/4,-3/8)) + 0(13~5) 

12*13 + 4*13~2 + 10*13^3 + 9*13-4 + 0(13-5) 

8. Future work 
Below we discuss some natural questions arising from our work. 

8.1. Global height pairings. Ultimately, we would like to compute the global 
height pairing. To do so, we would again require C to be a curve over a number 
field K with good reduction at each place v dividing p. We would also need a 
continuous idele class character \ '■ A^/if* — > Q p and a splitting H X (C '/ 'K v ) — 
H 1 ' ^/ K v ) © W v for each place v dividing p. Computing the local heights at those 
primes v away from p and those above p, the global height would then be the sum of 
all local heights. When K = Q, the recent Ph.D. thesis of Miiller [Mii llO addresses 
these local heights away from p, and putting together our results, we are able to 
compute global heights, as shown in Section [7J It would be quite interesting to 
extend these computations to number fields. 

8.2. Optimizations. In another direction, it is also of interest to optimize the 
present algorithm. Currently, the most expensive step is in computing the Laurent 
series expansion of a in the various Weierstrass local coordinates to reasonably high 
precision. As we are just interested in the residue of a J (3, is there a way to make 
this faster? 

8.3. Comparison with the work of Mazur-Stein-Tate. When the curve is 
elliptic, we are able to compare our algorithm for the height pairing with the al- 
gorithm of [MST06) , as in Section [71 We note that we compute the height pairing 
for divisors with disjoint support. It is obviously possible to compute without this 
assumption by replacing one divisor by a linearly equivalent one with this property. 
But it is also possible to extend the method described in |Gro86[ §5] . This extended 
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method can be compared directly with the method of [MST06] , as the height is just 
the height pairing of a divisor with itself. 

References 

[BBK10] J. S. Balakrishnan, R. W. Bradshaw, and K. S. Kedlaya, Explicit Coleman integration for 
hyperelliptic curves, Algorithmic Number Theory (G. Hanrot, F. Morain, and E. Thome, 
eds.), Lecture Notes in Computer Science, vol. 6197, Springer, 2010, pp. 16-31. 

[BcsOO] A. Besser, Syntomic regulators and p-adic integration. II. K2 of curves, Proceedings 
of the Conference on p-adic Aspects of the Theory of Automorphic Representations 
(Jerusalem, 1998), vol. 120, 2000, pp. 335-359. 

[Bes02] , Coleman integration using the Tannakian formalism, Math. Ann. 322 (2002), 

no. 1, 19-48. 

[Bes04] A. Besser, The p-adic height pairings of Coleman-Gross and of Nekovdf, Number The- 
ory, CRM Proceedings & Lecture Notes, vol. 36, American Mathematical Society, 2004, 
pp. 13-25. 

[Bes07] A. Besser, On the computation of p-adic height pairings on Jacobians of hyperelliptic 

curves, Sage Days 5, http://wiki.sagemath.org/days5/sched, 2007. 

[BeslO] , Heidelberg lectures on Coleman integration, Preprint, 2010. 

[CdS88] R. Coleman and E. de Shalit, p-adic regulators on curves and special values of p-adic 

L-functions, Invent. Math. 93 (1988), no. 2, 239-266. 
[CG89] R. F. Coleman and B. H. Gross, p-adic heights on curves, Algebraic Number Theory - in 

honor of K. Iwasawa, Advanced Studies in Pure Mathematics, vol. 17, 1989, pp. 73-81. 
[Col82] R. F. Coleman, Dilogarithms, regulators and p-adic L-functions, Invent. Math. 69 

(1982), no. 2, 171-208. 

[Col85] , Torsion points on curves and p-adic abelian integrals, Ann. of Math. (2) 121 

(1985), no. 1, 111-168. 

[Col89] , Reciprocity laws on curves, Compositio Math. 72 (1989), no. 2, 205-235. 

[Gro86] B. H. Gross, Local heights on curves, Arithmetic geometry (Storrs, Conn., 1984), 

Springer, New York, 1986, pp. 327-339. 
[Gut06] I. Gutnik, Kedlaya's algorithm and Coleman integration, M.Sc. Thesis, 2006. 
[HarOS] D. Harvey, Efficient computation of p-adic heights, LMS J. Comput. Math. 11 (2008), 

40-59. 

[Kau99] I. Kausz, A discriminant and an upper bound for u> 2 for hyperelliptic arithmetic sur- 
faces, Compositio Math. 115 (1999), no. 1, 37-69. 

[KedOl] K. S. Kedlaya, Counting points on hyperelliptic curves using Monsky-Washnitzer coho- 
mology, J. Ramanujan Math. Soc. 16 (2001), 323-338, erratum ibid. 18 (2003), 417-418. 

[Kob98] N. Koblitz, Algebraic aspects of cryptography, Algorithms and Computation in Mathe- 
matics, vol. 3, Springer- Verlag, Berlin, 1998, With an appendix by Alfred J. Menezes, 
Yi-Hong Wu and Robert J. Zuccherato. 

[MST06] B. Mazur, W. Stein, and J. Tate, Computation of p-adic heights and log convergence, 
Doc. Math. (2006), no. Extra Vol., 577-614 (electronic). 

[MT83] B. Mazur and J. Tate, Canonical height pairings via biextensions, Arithmetic and geom- 
etry, Vol. I, Progr. Math., vol. 35, Birkhauser Boston, Boston, MA, 1983, pp. 195-237. 

[MTT86] B. Mazur, J. Tate, and J. Teitelbaum, On p-adic analogues of the conjectures of Birch 
and Swinnerton-Dyer, Invent. Math. 84 (1986), no. 1. 

[MiillO] J. S. Miiller, Canonical heights on Jacobians, Bayreuth Ph.D. thesis (2010). 

[Nek93] J. Nekovaf , On p-adic height pairings, Seminaire de Theorie des Nombres, Paris, 1990- 
91, Birkhauser Boston, Boston, MA, 1993, pp. 127-202. 

[S+10] W. A. Stein et al., Sage Mathematics Software (Version 4-5), The Sage Development 
Team, 2010, http://www.sagemath.org. 

[Sch82] P. Schneider, p-adic height pairings I, Invent. Math. 69 (1982), no. 3, 401-409. 

Department of Mathematics, Massachusetts Institute of Technology, 77 Massachusetts 
Avenue, Cambridge, MA 02139 

Department of Mathematics, Ben-Gurion University of the Negev, P.O.B. 653, Be'er- 
Sheva 84105, Israel 



